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1. Introduction 


This project is the result of a joint effort between the 
Department of Electrical Engineering at the City Colleger CUNY; 
and the Department of Computer & Information Science at the 
Brooklyn College, CUNY. 

The goal of this project is to develop a software simulator to 
help NASA in the desiqn of the LMSS. The simulator will be used 
to study the characteristics and implementation requirements of 
the LMSS ' s configuration with specifications as outlined by NASA. 

This report represents progress made in implementing the 
simulator during the period 4/21/82 - 1/25/83. 
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2. The Present Simulator 


The basic components of the present simulator were described 
previously in the first report (NASA Grant: NAG 3-119). However, 
previously, not all the components of the system were operation'll 
and only scenario 1 was implemented. 

The current simulation now includes all the components, such 
as the interference sources, the S-band channels, the noise 
sources, the gateways, a revised master controller for the 
simulator, as well as other components which were recently 
completed . 

The present simulator is now capable of executing a.ll 5 
scenarios, which are described below . 

2.1 The Communication Scenarios 

Through the use of the five scenarios, it is posssible to 
study all the different possible combinations of types and modes 
of calls. The five scenarios are described below: 

(1) . Single Hop System: 

For those calls which require the communication between two mobile 
units, without going through a gateway, we group them together and 
form scenario 1. In this scenario, a mobile v/ 511. generate a call, 
for which the call will go up to the satellite, and be transponded 
down directly to the destination mobile. This can be seen in 
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figure 2.1. OF SW 

(2) . Double Hop, Single Gateway System: 

For those calls which require communication between two mobile 
units by going through a gateway, we group the together and form 
scenario 2. In this scenario, a mobile will initiate a call, 
which will be transponded by the satellite to a gateway. After 
processing by the gateway, the call is once again transponded by 
the satellite to the destination mobile. This can be seen in 
figure 2.2. 

(3) . Double Hop, Double Gateway System: 

For those calls which have to go through two gateways, we form 
scenario 3. In this scenario a mobile will generate a call, which 
is transponded to a gateway. The gateway towards the call to the 
master control station. The master control station then routes 
the call to the destination gateway. Then the call is once again 
transponded by the satellite to the destination mobile. This can 
be seen in figure 2.3. 

(4) . Mobile-to-Wirel ine System 

For any call that was initiated by a mobile and was desired to go 
through a fixed phone or to a mobile serviced by the cellular 
system, we have scenario 4. In this scenario, the call is 
initiated by the mobile and transponded by the satellite to a 
gateway. Then, the gateway towards the call into the wireline 
This can be seen in figure 2.4. 


network . 
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(5). Wire.line-to-Mobile System: 

If a call is iniated by a fixed phone or mobile in a cellular 
system t< a mobile in the LMSS,we resort to scenario 5. In this 
scenario a call is received from the wireline telephone network 
into a gateway. The call is then transponded by the satellite to 
the destination mobile. This can be seen in figure 2.5. 
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2.2 Output Facilities 

Also new in the present simulator is the revised version of 
the Master Controller (see 4.1.2, first report). In the present 
simulator it is possible to go through a range of values for any 
one of the components in the scenario. It is possible to consider 
a range of values for up to 14 different components 
simultaneously. Previously, if a type 2 measurement was chosen, 
the simulator would display a two column table of the input vs 
output SNR values, for a given set of system parameters. With the 
present simulator it is possible in a type 2 measurement for the 
simulator to display up to 15 columns (15 th column is output 
signal to noise ratio) , one for each range of component values. 
In the present simulator, it is also possible to get plots of the 
input vs. output SNR values for any range of values of the other 
components (see section 5, "Sample Runs"). The looping or ranging 
feature mentioned above can also be used in a type 1 measurement, 
were the input signal, output signal (no disturbances present), 
and output signal (with all specified disturbances) are displayed 
in a table. It is now possible to consider a type 1 measurement 
over a range of values (up to 14 possible ranges of parameters) . 
Once again, a plotting facility is available for every type 1 
measurement table displayed (see section 5, "Sample Runs"). 

The revised master controller function has the following 


algorithm : 
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******************************************************************* 

I 

/SIMULATION MASTER CONTROLLER FUNCTION/: j 

initialize simulation master controller function> ! 

DO <range of Specular-to Multipath Power Ratio — UHF Uplink > 

l 

<range of Specular-to-Multipath Power Ratio — S-Band Downlink> i 

<range of Specular-to-Multipath Power Ratio — S-Band Uplink^ 

<range of Specular-to-Multipath Power Ratio — UHF r-ownlink> »j 

i 

<range of carr .-to-interf . power ratio — UHF Uplink Interference> j 

I 

<range of phase values--UHF Uplink Interference *> = 

<ranqe of carr , -to-interf . power ratio — S-Band Downlink Interf> j 

i 

<range of phase values — S-Band Downlink Interference> 1 

<range of carr .-to-interf . power ratio — S-Band Uplink Interf.> 

<range of phase values — S-Band Unlink Interf.* 

<range of carr . -to-interf . power ratio — UHF Downlink Interf> 

<range of phase values— UHF Downlink Interf.* 

<RANGE OF SNR > 

<calculate SNR> 

<calculate standard deviation for white noise> 

<perform simulation reset> 

<set flag for mobile transmitter> 

<PERFORM a simulation run> 

<perform master controller instrumentation package> 

END <signal to noise ratio> 

END /simulation master controller function/ 
******************************************************************* 
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2.3 System Structure 

The present modular structure can be seen in the flowchart of 
figure 2.6, with its subroutines defined in table 2.7. 



F igure 2.1 ; Single Hop Single Gateway System 
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3. Processing Voice Via The Simulation 


One of the systems greatest present attributes is its ability 
to process voice. The present simulator can be utilized to listen 
to any of the following effects on voice (via each of the five 
scenarios) : 


(a) . 

Listen 

to 

the 

effects 

of 

noise . 

(b) . 

Listen 

to 

the 

effects 

of 

interference. 

(c) . 

Listen 

to 

the 

effects 

of 

modulation. 

(d) . 

Listen 

to 

j-i 

UIC 

« jl i. c:C uo 

o f 

the fading channel. 

(e) . 

Listen 

to 

the 

effects 

o f 

filtering on the voice. 

(f) . 

Listen 

to 

the 

effects 

of 

the space diversity receiver. 

(g) . 

Listen 

to 

the 

effects 

of 

the satellite nonlinearity. 

(h) . 

Listen 

to 

the 

effects 

of 

varying sampling rates. 

(i) . 

All combinations of the . 

above . 


An analysis of a number of the above tests are included within 
the report (see section 6.2). 

A basic model of how to use the simulation to process voice is 
shown in figure 3.1, In this model, an analog to diqital 
converter samples voice at a rate of at least 8 kilohertz (minimum 
sampling rate for voice bandlimited to 4000 Hz.), producing 12 bit 
samples. In the second stage, these 12 bit samples are converted 
into 16 bit integer format (IBM format) . These 16 bit integers 
are then read by the LMSS simulator. (see section users manual 
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for simulator voice setup) . The output of the LMSS simulator is 
then converted from 32 bit floating point form (IBA format) back 
into a 12 bit integer format by the fourth stage. In the fifth 
stage, a digital to analog converter with the proper hardware 
connected , converts the samples into voice. In between each of 
the above stages, any type of stornqe device could have been used 
for intermediate storage. 

The system outlined above is general, and actual 
implementation depends on the particular intermediate systems 
involved (stages 1,2, and 5). 

In figure 3.2 we have a detailed block diagram of the specific 
algorithim we used in implementing what was shown in figure 3.1. 

The first stage consisted of a 12 bit analog to digital 
converter interfaced into a Hewlett Packard 1000 (HP 1000) 
computer system. This system sampled the voice (band! imited to 
4000 Hz) with a sampling rate of 0 Khz, and stored it on tape with 
a "WRITT" format (This is an HP utility which stores data onto a 
tape in an ultra-condensed coded format, which can only be read by 
another HP system) . We had to store it in this form because the 
system we used didn't have enouqh storage focilty to do it any 
other way. In the next stage, we used another HP 1000, with much 
greater storage facility to convert the WRITT format tape into 
straight ASCII data. For error control and checking purposes, 
data was stored in the format of an integer (1R bit integer, 
represented by 4 characters) followed by a space. 

This tape was now to be mounted on the IBM machine. We used 
the program shown in figure 3.3 to convert from ASCII into EBCDIC 


OF PCOJ5 QUALITY 


(IBM) format. We stored the data on a disc. We than used the 
LMSS simulator to process the voice samples. 

The LMSS simulator has the facility to read in integers 
(16-bit) in the format mentioned above. The simulator then 
interpolates between these samples if specified by the user, and 
then processes them. The output data set is a stream of 32-bit 
(IBM format) floating point numbers. Module 3.4 (see figure 3.4) 
is then used to convert this 32-bit data stream into a 16 bit 
interger data stream (format 16(lx,I4)), with sample values 
normalized to plus or minus 2047 (12—bit number). Module 3.S (see 
figure 3.5) is then used to convert the data from EBCDIC to ASCII 
format, and stores it on tape. The tape is then read in by module 
3.6 (see figure 3.6) using the HP-1000. This module converts data 
into the format needed by the digital to analog converter. Then 
finally the D/A converts the samples back to voice. 

All the above steps were needed because of code conversions 
between machines as well as numerical format (the IBM is a 32-bit 
machine while the HP-1000 is a 16-bit machine) . Ideally, one 
would use the setup of figure 3.1 if the proper equipment is 


attainable 


VOICE 






























of spoon (WAslsbv 


Thi$ proaraw reaas n tape m ASCII Format, arm 
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//JOB 

//STEP! EXEC PGrr-PREU!* 

//SYSPRINT I)D SYSOIJT - A 
//HYB IN DI> DUMMY 

//SYSIJT 1 DI) DBM-Owr , D T BP- ( Oi D PARS ) , LABEL" ( j. , Ml. ) , 
// I.INi 'l - I APE ,DUP.-ll KICI -HO , p.l. KB ) '/h - HO , Kr ( :| M- (■ R , J)| 
//BYSIJT2 DD I)BM = UiYI. . CIC „ BMC . TAB I? I) , UN] I I O , I) J BP= < 
// VUI. “Si: RaftCKOO? , DUB- ( I I?) Cl -H2MM , P.l KB j ;■(- ^ K23B , |y(. 
/ / S PACE = ( trk , ( 20 , fi ) , Ml BE" ) 

/>< 

// 


o o n u t* v % s 3 t t o 


VUL ~BER"CCT2BO , 
W«M , UP I CD-U ) 
rUATI.R) , 

I.Hi-ll , DBURG-PJ5 > 


(27 


o^u'f e 


5 .*+ 
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rr.* 

** fcCjJ 


-1- 

4. 

5. 

6 . 

7. 

8 . 
8 . 

40, 
11 . 


+ t,i> MJ.HWI ..,*4 r-tfiine, a Fiji* at-' 32 bit <•' 1 KWt- J. tf*t i'HUtt . HMMU.OfS 

n ijviiitui ' k(ii/l,y> i»«r‘ ri!i:ori 0 . aim oomjf, '/?*.*» tnofii xvito a * uc m - 1** 
HIT. 1 1 1 ( (11 D 0 T' *3 ( 3 Vi t tt Si 01* { » ), F lit’ UM» in tilt* D/A *5 Y ** t P til • 


// JOB TIME«4,RERION“400K 
//•ttMAIN LINES = 50 
// EXEC FORTGCQ 

//FORT-, 8YB1N 1)1) * ' ' 

C •«. it % » •)(• ### -It ^ It IHHt It It «■ 41- •), It II -It -It -)(- * It It -It It It It II II It » -it- 4HI # It # * ## ■# «• » It It 'll It It It It It It It It It H- IHHt It 


12. 


Rtf At MAX 

13. 


r.) I MENS I OH X ( 82768 ) , I X ( 32768 ) 

14. 


MAX = 0 

15. 


HEAD (1,1 00 , END-200 ) ( X ( J ) , J " 1 , 32 768 ) 

-1G. - 

3 00 

FORMAT <4 (If 14.7) ) 

17 p 


DO 55 J® 1,3 2 768 

18. 


l p ( ABS ( X ( J ) ) . GT . MAX ) MAX = ABI5 ( X < J ) ) 

19 . 

l- « 

COIvilHMI- 

-20 « 


PRINT 70, MAX 

21 

70 

FORMAT ( ' ' , ' MAX® ' ,E14.7) 

-Z-2-* - 

... ... . . 

DO 66 J= 1r3Z7G0 • . • 

23. 


X ( J ) =X ( J ) /MAX 

24 . 


X ( J ) - X < J ) It 2 046. 

25. 


X< J> «X( j >+2u*7. 

-26*- ■ 


X ( J ) -X <.))+, 5 

27. 


I X ( J ) =>(< J ) 

-28. 

KB 

CO NT I HUE 

29. 


OR 116(7, 333 ) ( I X ( J ) , J - 1 , 32768 ) 

30. 

\ j t!»»7 

F ORMAT ( 1 0 ( 1 X , T 4 ) ) 

31 . 

zoo 

CONTINUE 

w 


STOP 

33 • 


END 

-34 „ 

//GO. 

PTOl FOO 1 DI) DSN-WYL . CC . BHC . RB1B , D.ISP-SHR 

35. 

//GO. 

F | 071- 00 1 1)1) DHiv:- WYL . CC . PUB . XS i. B , \M I f = 3330 , 

36. 

// D ISP* < KlEW r CAT l G r CAT LG ) r RPACF- (TRKr(14r5>rRI wb ) r 

37 . 

// 

001. « BE R= SCR 003 , 

3B . 

// 

1)CB“ ( 1)HURG-*PS , Rt-CFN-FB , L REEL .“80 , 81. KR I 76 = 6,1 BO ) 

39. 

//GO. 

SYS IN 1)1) it 


tenure. . 
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// JOB T1ME = S r Rf" G I ON-BOOK 
//8TEP1 EXEC HI in u lb BUI- Mbit 
//bYBFRlWl 1)1) SYSlJlIT^A 
//SYS IN DD DUMMY 

//SYSUT1 DD DSM-UYI . CC . PUB . F Il.Kl 1 » D 1. SP=SHR t 
// I IN i i -HYHDA 

//SYSUT2 DD DSN^ONt ,tlNI i « TAPE , DCB* ( LRECL «80 r Bl.KS I ZE = 80 r RECF'M«FB r 
// DPN=B» (IP i f'D-N ) i Mill -Shrt-n iOO'CI AIM l -< 1 , IMI ), D 1 SP = < NEW r KEEP ) 

/•«• 

//STEPS EXEC PGM- PRESS 
//SYS PRINT DD SYBOIJT-A 
//SYS IN DI) DUMMY 

//SYSI I I 1 Dl) DSKo- flNE ,i)l SP» ( OLD , KEEP > » LABEL* < 1 rNl. ) , V0L«BER*CCT004 , 
// UNI I -TAPP r DCS- i I RE l!l -SOrBI KS I 7\ -HO , Rl- Cl- M=EB r DEN-3 , 0PTCD = (3 ) 
//SYSUT2 I)D DSN- NYT . CC . PUB . NASA , UN Li *SSSO , D) SP- ( , CATl G ) , 

/ f U 01“ S b R - S C R 0 O 2 i- D ( : B - ( L R E C I -SPSS , Si K S i 7E-G7SS r RECF- I r DSORG - PS ) r 
/ / S PAl.ir — ( 'I RK , ( >\ n<> , S ) , Kl SH ) 
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4. User's Guide 

This user's guide will provide the user with a step by step 
explanation of the questions that will be asked during a typical 
session with the simulator. References are given in each question 
so that the user can refer to the appropriate sections if an in 
depth discussion about the topic is desired (references made refer 
to the first report, "A Computer Simulator for a Mobile Telephone 
System", NASA Grant 3-119). 

The following is a typical question-and-answer session with 
all the possible questions either being asked or presented in the 
explanation. The answers shown here are for a batch version 
version for a run. In an actual run, the answers shown would come 
out on the same line as the question. However we have seperated 
them here onto seperate lines for clarity. For an online session, 
merely the question will be prompted unto the user terminal, for 
which he would interactively enter his response. 

The first question asked prompts the user as to whether the 
session is an online or batch session, and is displayed as 
follows : 

Is this an online session? (Y/N) (1) 

N 

The user should answer a 'Y' if it is an online session, or an 
'N' if it is a batch session. In this case the response is "N" 
(no^ . 

Next, the type of calls that are available in the 
communication system are displayed as follows: 
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communication system are displayed as follows: 

TYPE OP CALL AVAILABLE: (2) 

1: M1->M2, rural mobile to rural mobile in same UHP beam 
2: M1->M3, rural mobile to rural mobile in different UHF beam, 
in same S-band 

3: Ml->M4, rural mobile to rural mobile in different UHP beam, 
in different S-band. 

4: M1->FC1 , rural mobile to fixed in same S-band beam 

5: M1->FC2, rural mobile to fixed in different S-band beam 

S: FC1->M1, fixed to rural mobile in same S-band beam 

7: FC2->M1, fixed to rural mobile in different S-band beam 

INPUT TYPE OF CALL TO BE SIMULATED: 

1 

Here the type of calls that are available in the 
communication system are displayed, each one being designated by 
a number from 1 through 7. Right after the list, the user is 
prompted to input the type of call he wishes to be simulated by 
inputing a number from 1 through 7. In this case, the first call 
("l”) was chosen. For an explanation of each of these calls, see 
section 2.2. 

After the user has inputed the type of call to be simulated, 
the modes of call available within the type of call selected will 
be displayed (refer to section 2.3). In this case, the user chose 
to simulate type 1, which is rural mobile calling rural mobile in 
the same UHF-beam and S-beam. Correspondingly, the following 
choice of modes was displayed: 

MODE OF CALL AVAILABLE: (3) 

1: M1->M2 , hard wired transponder 

2: M1->M2, direct switched transponder 

3: Ml->M2, indirect switched transponder 

4: M1->G1->M2, double hop system 

INPUT MODE OF CALL TO BE SIMULATED: 

1 

Note that the type of modes that will be displayed depends on 
the type of call that is chosen (in this case type 1 was chosen). 
A full list of the range of possible modes corresponding to the 
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type of call chosen is found in section 2.3. (table 2.3.1). 

One notes that following the list of modes available, the user 
is prompted to input the mode of the call by entering its 
corresponding number. In this case the user chose w l", from a 
choice of l to 4. 

Based on the type of call an. mode of call selected, the 
program will determine which scenario to be simulated (refer to 
section 3.3). The questions asked from this point on will vary 
depending on the features to be simulated within the scenario. 
Even though each of the five scenarios are different, the 
questions to be answered for each of them are more or less of the 
same nature. Therefore the rest of our discussion will focus only 
on those questions that are asked when the single hop system is to 
be simulated (which is the case of this par.icular simulation 
run) . 

In the next question, the user is asked to choose from the 
types of signal sources available. Choice 1, corresponds to a 
program generated sinusoid which can be used for making analytical 
tests. If choice "1" is chosen, no further questions aout the 
signal source are asked and question (4) is then displayed. 
Choice "2” corresponds to sampled voice. If type 2 is chosen, 
then see section 3.2 for an explanation of procedures to be used. 

TYPES OP SIGNAL SOURCES AVAILABLE: 

1: Program generated single tone sinusoid. 

2: Sampled voice from tape source. 

Choose type of signal to be used (II): 

1 

In the next question, 

FREQUENCY OP THE BASEBAND SIGNAL ( LESS THAN 3000 HZ . 1 (4) 

(IN HERTZ, F7. 2) : 


1000.00 HZ 
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U ? Te user is asked to input the frequency of the baseband signal in 

hertz using the format "F7.2". The baseband signal has an upper 

limit of 3000 Hz (refer to section 5). 

The next question displayed is as follows: 

NOTE: POWER OF THE BASEBAND SIGNAL HAS TO BE LESS OR BOUAL (5) 
THAN. 5 WATTS; 

SO THAT THE MAXIMUM INSTANTANEOUS FREQUENCY DEVIATION IS 
LESS OR EQUAL THAN 12000 HZ 

POWER OF THE BASEBAND SIGNAL (IN WATTS, F7.2): («) 

0.50 WATTS 

Here the user inputs the power of the baseband signal in watts 
using the format F7.2. The power of the baseband signal has to be 
less than or equal to .5 watts, so that the maximum instantaneous 
frequency deviation is less than or equal to 12000 Hz (refer to 
section 5) . 

In the next question that appears 

THE CARRIER POWER (IN WATTS, F7.2): (7) 

1.00 WATTS 

the user is asl<$d to input the carrier signal power in watts using 
the format F7.2. 

The next question appearing, 

THE FREQUENCY DEVIATION (IN HERTZ, F8.2): 12000 HZ. (R) 

(FIXED FOR NOW) 

asks the user to input the frequency deviation in Hertz using the 

format ,F8.2. Under the current version, the frequency deviation 

* a fixed at 12,000 Hz in order for the simulator to meet AMPS 

specifications (refer to section 5). 

The next question is displayed as follows: 

How many times the nyquist rate do you want the sampling (9) 
frequency to be" 5 (2 - 4, II): 
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Here the user is asked to input fchs campling frequency in terms of 
a multiple of the Nyquist rate. Acceptable values are from twice 
the nyquist rate to four times the Nyquist rate. Twice the 
Nyquist rate is the minimum in or er for the simulator to 
intergrate and differentiate correctly. Four times the Nyquist 
rate is the maximum since any increase of the samplinq frequency 
beyond this point will not improve any significant approximation, 
but will prolong run time substantially (refer to section 3.2). 
tn the next question, 

THE CARRIER FREQUENCY (IN HERTZ, F7.?):(N0T USED FOR NOW) (10) 

the user is asked to input the carrier frequencyin Hertz using the 

format F7.2. Under the current version of the simulator, the 

carrier frequency is not used in the simulation. 

The next question is displayed as follows: 

DO YOU WANT THE COMPRESSGR/EXPANDER IN? Y-YSS , N-NO : (11) 

Y 

Here the user is asked whether the compressor/expander (AMPS 
specs) are used in the transmitter and receiver (refer to section 
5) , This question is answered either with a "Y" (yes) or an "N" 
(no) . 

The next question appears as follows: 

DO YOU WANT THE PRE-EMPHASIS - DE-EMPHASIS FILTERS IN ? (12) 

Y-YES , N-NO : 

Y 

Here the user is asked whether pre-emphasis and de-emphasis 
filters (AMPS specs) are used in the transmitter and receiver (see 
section 5 ) . 

The user is then prompted as to whether or not interference is 
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present in the UHP uplink sattelite channel. This question is 
displayed as follows: 


DO YOU WANT A CO-CHANNEL INTERFERER IN THE UPLINK? 


Y-YES,N-NO: 

Y 


Are multiple values of the carrier to interference power 
ratio to be tested? (Y/N) : 


Y 

Input the range of values to be tested (in DB) : 
Input the initial value (-99 — +99,13): 

10 DB 

Input the increment value (01 — 99,13): 

10 DB 


Input the final value (01 — 99,13): 

70 DB 


Are multiple values of the interference phase values 
to be tested? (Y/N): 

Y 


Input the range of values to be 
Input the initial value (FH.2) : 

Input the increment value (F6.2) 

Input the final value (FH.2) : 


tested fin degrees): 
10.00 DEGREES 
5.00 DEGREES 
25.00 DEGREES 


( 13 ) 

(14) 

(15) 

( 16 ) 

(17) 

( 18 ) 

(19) 

( 20 ) 
( 21 ) 


Note the answer to the first of these questions was "Y" (yes) . 
If the response was "N" (no interference) , then the next eight 
questions would not have appeared. However, since in this case 
the response is yes, they did appear. In the first of these 
questions (14), the user is asked if he wishes to input a range of 
values for the carrier to interference power ratio. If the answer 
is "Y" (yes), as in this case, then questions 15-17 will follow. 
Here, the user inputs the initial value, increment value, and 
final value respectively (in dB, using the format 13) of the 
carrier to interference power ratio. If the answer to question 
(14) was "N" (no), then the user would have just been asked to 
input one value of the carrier to interference power ratio in the 
same format as above. After question 1?, the user is asked wether 


he wants to 
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input a 

range 

of 

phase 

values 

for 

the interference source 

(question 

18) . 

If 

the 

response 

to 

this is "Y H 

(yes) , then 

questions 

19-21 

will 

appear 

, asking 

the 

user to input 

the initial , 


incremental, and final value of the phase, respectively (using 
format Ffi»2 , in degrees). If the response to question 18 was H N H 


(no), then a question would have appeared aakinq only one value 
for the phase. 

In the next question, 

IS FADING PRESENT IN THE UPLINK CHANNEL? Y-YES,N-NO: (22) 

Y 

one responds as to whether or not the uplink satellite channel is 
a Fading channel or not. One responds to the question by either 
typing a tt Y" for yes or a "N" for no. If the response is no, the 
next five questions (questions 28-28) will not appear. However, 
as in this case, if the response is yes, then questions 24 through 
28 will definitely appear. If interference is present, as Is the 
case in the current simulation run, then question 23 will also 
appear. If interference were not present in the uplink, then 
question 23 would have not been asked. Questions 2,3-28 are 
displayed as follows: 

IS THE INTERFERENCE FADED? Y-YES,N-NO: (23) 

Y 

THE UPLINK FADING CHANNEL IS PRESENT 

TYPES OF FADING CHANNELS AVAILABLE: (24) 

1: NO SPECULAR COMPONENT (RALEIGH FADING) 

2: SPECULAR COMPONENT, SHORTEST PATH 
3: SPECULAR COMPONENT, MEAN PATH 
INPUT TYPE OF CHANNEL: 

2 

ENTER THE MULTIPATH SPREAD TIME (IN MICROSECONDS, F9.2): (3S) 

sno.no MICROSECONDS 

ENTER THE DOPPLER SPREAD BANDWIDTH (IN HERTE , F7.2): PR) 

1.00 HR 
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Are multiple values of the specular to multipath power (27) 
ratio to be tested? (Y/N) s 

N 

ENTER THE SPECULAR-TO-MULTIPATH POWER RATIO (IN DB , 13): (28) 

100 DB 

In the first of these questions, the user respond's as to 
whether or not the interference is effected by the fading channel. 
If it is, it is done so independently of the signal and is added 
at the end. If it is not, then the interference is just added to 
the faded signal. 

In the second of these questions, the user decides which type 
of fading is present, by entering a 1 digit number from 1 to .3 
corresponding to the type of fading that is present. If (1) is 
chosen, then the fading channel will be a Rayleigh distributed one 
with no specular component. If (3) is chosen, then the channel 
will have Rician statistics, with the specular component in the 
middle. In other words, the mean path is taken. In the last 
case, (2), the channel has Rician statistics with the shortest 
path being taken (refer to section 6) . 

In the third of these questions (25), the user is asked to 
input the total multipath spread time of the fading channel in 
microseconds using the format F9.2. In this case, "500.00" was 
entered as data. 

In question (25) , the user is asked to input the doppler 
spread bandwidth of the fading channel in Hertz using the format 
F7.2. In this case the number "1.00" was entered as data. 

Question (27) will only be asked if the response to question 
(24) is either 2 or .3. If the response is 1, then Rayleigh fading 
is being used which exhibits no specular component. Question (2?) 
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once again asks wether or not the user wishes to input a range of 
values for the specular to multipath power ratios (in dB) . In 
this case the answer is "N" (no) , and thus question (28) appears , 
which asks the user to input only one value of the specular to 
multipath power ratio (in dB) using the format 13 (in this case 
100 was entered as data) . If the response to question (27) was 
yes, then three questions asking for the initial, incremental, and 
final values respectively (using the same input format as question 
(28)) would have been displayed as they have in the questions 
asked previously. 

The next question. 


DO YOU WANT A CO-CHANNEL INTERFERER IN THE DOWNLINK' 5 (29) 

Y-YER , N-NO : 

Y 

Are multiple values of the carrier to interference power (30) 
ratio to be tested ? (Y/N) : 

N 

INPUT CARRIER TO INTERFERENCE POWER RATIO IN DB(I3): (31) 

30 DB 

A-e multiple values of the interference phase values (32) 

to be tested ? (Y/N): 

N 

INPUT PHASE OF INTERFERENCE IN DEGREES (F7.2): (32.1) 


30.80 DEGREES 

is the dual of the uplink interference case. The data entry is 
exactly the same as it were for the uplink case. Note however, in 
this case that the responses to the questions asking if the user 
wants to input a range of values are both negative. In this case 
the user is only inputing single values (030 dB for the power 
ratio, and 30.80 Degrees for the phase). 

The next question, 

IS FADING PRESENT IN THE DOWNLINK CHANNEL? Y-YES, N-NO: (33) 

Y 


t 


is basically the dual of the uplink case, the only exception 
however is the facility for a space diversity receiver. Once 
again, if the response to the above question is "N" (no) , then 
questions 34 through 45 will not be asked. If the response is "Y" 
(yes), the next question that will be displayed is: 

IS A SPACE DIVERSITY RECEIVER BEING USED? Y-YES,N-NO: (34) 

Y 

If a space diversity receiver is being used (response "Y") , as 
is the case for this example, then parameters for 2 fading 
channels will have to be set; questions 35 through 40 for the 
first channel and questions 41 through 45 for the second channel. 
If a space diversity receiver is not being used (response "N") , 
then the user will only be prompted to enter parameters for one 
channel . 

If co-channel interference is present in the downlink channel 
(response to (29) is "Y") , then the following question will be 
displayed : 

IS THE INTERFERENCE FADED? Y-YES,N-NO: (35) 

Y 

Here one enters whether or not the interference should be 
faded. If yos (as in this case), the user enters a "Y", and if 
not, the user enters an "N" . 

Questions 35 through 40 or questions 35 through 45 if a space 
diversity receiver is present, are entered in an analagous way as 
were the questions for the uplink satellite channel. 

THE DOWNLINK FADING CHANNEL IS PRESENT 

SET PARAMETERS FOR FIRST FADING CHANNEL: (35) 

TYPES OF FADING CHANNELS AVAILABLE: 

1: NO SPECULAR COMPONENT (RALEIGH FADING) 

2: SPECULAR COMPONENT, SHORTEST PATH 
3: SPECULAR COMPONENT, MEAN PATH 
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INPUT TYPE OP CHANNEL: 

2 

ENTER THE MULTIPATH SPREAD TIME (IN MICROSECONDS , F9.2): (37) 

500. 00 MICROSECONDS 

ENTER THE DOPPLER SPREAD BANDWIDTH (IN HERTZ, F7.2): (38) 

1.00 HZ 

Are multiple values of the specular to multipath power (89) 

ratio to be tested? (Y/N) : 

N 

ENTER THE SPECULAR-TO-MULTIPATH POWER RATIO (IN DB , 13): (40) 

100 DB 

SET PARAMETERS FOR THE SECOND FADING CHANNEL: (41) 


TYPES OF FADING CHANNELS AVAILABLE: 

Is NO SPECULAR COMPONENT (RALEIGH FADING) 

2: SPECULAR COMPONENT, SHORTEST PATH 
3: SPECULAR COMPONENT, MEAN PATH 
INPUT TYPE OF CHANNEL: 

2 

ENTER THE MULTIPATH SPREAD TIME (IN MICROSECONDS, F9.2): (42) 

500.00 MICROSECONDS 

ENTER THE DOPPLER SPREAD BANDWIDTH (IN HERTZ, F7.2): (43) 

1.00 HZ 

Are multiple values of the specular to multipath power (44) 

ratio to be tested? (Y/N) : 

N 

ENTER THE SPECULAR-TO-MULTIPATH POWER RATIO (IN DB , 13): (45) 

100 DB 

In the next question. 

Are multiple SNR values to be tested? Y-yes,N-no: (46) 

Y 

the user is asked whether multiple SNR values are to be tested. 
If the response is "Y" (yes), then questions 4-7.1 through 4*7.3 
will appear asking the user to input the range of SNR's (in this 
case, the carrier to noise ratio) for the Gaussian noise in the 
mobile receiver. If multiple SNR values are not to be tested 
(response "N") , then the user will be asked to input only one SNR 
value . 

Questions 47.1/47.? would appear as follows: 
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Input the range of SNR values to be tested (In Db) : (47) 

Input the initial value of SNR (-99 — +99, 13): (47.1) 

36 DB 

Input the increment value (01 — 99, 12): (47.2) 

4 DB 

Input the ending value of SNR (-99 — +99, 13): (47.3) 


40 DB 


In the next question, 

IS A SPACE DIVERSITY RECEIVER REQUIRED? (Y/N) (4B) 

Y 

the user is asked whether or not a space diversity receiver is to 
be used. If it is present (response "Y") , then the user will be 

asked (next question) to input the duration of the decision period 
which is the duration of the lapse time before a decision is made 
to determine which receiver is receiving a stronger signal . If it 
is not present (response "N") , then the next question (49) will 
not displayed. 


INPUT DURATION BETWEEN DECISION TIMES FOR THE S.D. RECEIVER 
(F7.5, IN SECONDS) (49) 

0 . 00020SEC . 


The next question is displayed as follows: 

Input the approximate duration of simulation in seconds (50) 
(0.01 - 9.99, F4. 2) : 

0.01 SEC. 

Here the user inputs the approximate duration of the input signal 
to be tested in seconds using the format F7.5. For now, the 
boundary imposed is from 0.01 seconds to 9.99 seconds (refer to 
section 4.1). 

In the next question displayed, 

Type of performance measurement available: (51) 

1. Compare recovered output signal to original input signal 

2. Measure output siqnal to noise ratio 
Input type: 

1 
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the user is asked to input the type of performance to be measured. 

Under the current version, the user has the choice of choosing 

either to compare the recovered output signal to the original 
input signal (choice "1") , or (choice "2"} to measure the output 
signal to noise ratio vs input signal to noise ratio (refer to 
section 4.2). If choice 2 is chosen, the user will be asked 1 if he 
wishes that a plot be displayed of the input vs output SNR values, 
asked. However, if choice 1 is chosen, as is the case in this 
example, than the following will be the last question displayed: 

DO YOU WANT OUTPUT TO BE PLOTTED? (52) 

Y 

Here the user has the option of having the output vs input 

signals plotted. If the response to this question is "Y" (yes), 

then following the printing of the input vs output signal values 
(for each respective SNR value), a graph will be plotted 
corresponding to these values. The maximum possible length of the 
graph is 2000 lines (per SNR value). If the response is "N" (no), 
then just the values without a graph will be printed. 

A few lines of the resulting output are printed as follows: 
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Input signal Output signal (No Noise) Output signal (With Noise) 


The SNR is 

36. 

00 DB 




0 . 1000000E 

01 

-0 . 3008653E 

02 

0.6992087E 

02 

0.9972572E 

00 

-0.9639141E 

03 

0.2256/19 3 E 

04 

0.9890135E 

00 

-0.7375438E 

04 

0. 1755339E 

05 

0.9752213E 

00 

-0.2662992E 

05 

0.6575075E 

05 

0 . 9557958E 

00 

-0.6205 3 68E 

05 

0.1653216E 

06 

0 . 9306050E 

00 

-0.1081285E 

06 

0.3328434E 

06 

0 . 8994539E 

00 

-0.1487061E 

06 

0.5R98719E 

06 

0 . 8620588E 

00 

-0 . 1609989E 

06 

0.9725955E 

06 

0 . 8180044E 

00 

-0.1277967E 

06 

0. 1540994E 

07 

0.7666720E 

00 

-0.5651906E 

05 

0.23R7130E 

07 

0.7071075E 

00 

-0 .781881 3 E 

03 

0. 3643070E 

07 

0.6377603E 

00 

0.7991531E 

05 

0 . 5487024E 

07 

0 . 5558942E 

00 

0.4918778E 

06 

0.8145275E 

07 

0.4559756E 

00 

0.1515125E 

07 

0.U88825E 

08 

0.3 23311 2E 

00 

0. 3496630E 

07 

0. 1701875E 

08 

0, 1490450E- 

-02 

0.6825587E 

07 

0. 2384659E 

08 

-0. 3233029E 

00 

0 . 1189424E 

08 

0. 3265666E 

08 

-0. 4559686E 

00 

0 . 1904904E 

08 

0. 4366894E 

08 

-0 . 5558881E 

00 

0. 2853454E 

08 

0. 5699677E 

08 

-0 . 6377543E 

00 

0.4043526E 

08 

0.7258918E 

08 

-0.7071018E 

00 

0 . 546204 0E 

08 

0 . 9018307E 

08 


720. -0.81J4990E 00 -0.3250146E 05 -0.3252382E 05 
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5. Sample Runs 

One of the biggest advantages of the present simulator is its 
ability to span over a range of inputs for the different devices, 
as opposed to only entering one value at a time as was the case in 
the previous simulation. This facility is illustrated in figure 
1. In this particular example, 4 parameters are variables: the 
specular to multipath power ratio of the UHF downlink channel, the 
carrier to interference power ratio of the UHF downlink co-channel 
interfer, the phase values for the UHF downlink co-channel 
interferor, and the input signal to noise ratio of the Gaussian 
noise source. One can see from figure 1 that every possible 
permutation of the above values is exhibited. In this simulation 
run, a type one measurement is chosen with the plotting facility 
on. Thus the output is a list of the input vs. output values of 
the system, followed by its corresponding plot. One notes that a 
type one output exists for each of the possible permutations as 
mentioned above. 

Figure 2 illustrates the above mentioned "looping" feature for 
a type 2 measurement. In this output, a type 2 measurement of the 
input vs. output SNR values a id its corresponding plots are 
exhibited for every possible permutation of the 7- input 
parameters (columns) shown. 

In figure 3, the output SNR values are being plotted as a 
function of the input SNR values for different values of the 
carrier to interference power ratio of the UHF uplink interference 
source. Again, the looping facility is exhibited here. One notes 
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from the output curves, that as the carrier to interference power 
ratio gets smaller, the output curves flatten out. At large 
values, the curves remain linear above threshold as expected. 

Figure 4 illustrates the execution of a double hop system 
(scenario 2 in this case) . One notes here that both an S-band and 
an UHF band uplink channel are set, as well as an S-band and an 
(JHF band downlink channel. 
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The results of the present simulator are divided into two 
sections; The analytical results (section 6.1), and the processed 
voice results (section 6.2). In both cases, tests are made using 
a single hop system (scenario 1). 

6.1 Analytical Results 


The system discussed can be seen to have the responses shown 
on various plots obtained for several different Rician channels 
over a single hop system (scenario 1). 

We first analyzed a single hop system over a channel without 
fading or interference present as a function of the input SNR 
values. A pure sinusoid was used as input. The results can be 
seen in figure 6.1 (referenced from first report). Note that the 
results are analyzed for various combinations of signal processing 
staqes. These results represent a pure FM system with white noise 
present . 

We next simulated a single hop system over a Rician channel 
with no interference present, and analyzed the output vs input SNR 
values for various values of the specular-to-multipath component 
power ratios of the fading channel. The results are summarized in 
the plot of figure 6.2. First, we note that when the 
specular-to-multipath power ratio of the channel is 100 db, no 
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real effects of fading can be seen and the output vs input SNR 
curve remains a linear one as expected. An analysis of the 
various plots, however, does show that as the specular-to- 
mul tipath power ratio is decreased, a flattening out of the linear 
curve does occur and the effects of fading become increasingly 
destructive . 

We next analyzed a single hop system over a Rician channel 
with no interference present, with and without a space diversity 
receiver present . The results are summarized in figure 6.3. One 
observes from this plot that there is a significant improvement 
with the space diversiy receiver on. 

We next plotted the output signal-to-noise ratio as a function 
of the specular-to-multipath power ratio of the fading channel, 
with no other disturbances present. We started with a 
specular-to-multipath power ratio of -35 dB (practically Raleigh 
fading), and went all the way u,p to a specular-to-multipath power 
ratio of 100 dB. The results are summarized in figure 6.4. Here 
we noticed a threshold effect beginning at a specular-to-mul tipath 
power ratio of about +05 dB. The curve moves up very quickly in a 
linear fashion from 5 dB to 50 dB, and then flattens out to its 
maximum value. 

Figure 6.5 considers the effects of interference on the output 
signal to noise ratio. One can observe from the plot that an 
interference source with a carrier to interference power ratio 
ranging from 100 dB all the way down to 20 dB has no major effect 
on the output performance. However, as one gets below 20 dB, a 
flattening out of the input vs. output SNR curve does occur above 
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threshold. As the carrier to interference power ratio is taken 
below 8 dB, the flattening out effects become very predominant and 
output performance starts getting poor. At 0 dB, the curve 
flattens out completely to about .5 dB (Threshold effect). 

Figure 6.6 considers results obtained using scenario 2. In 
this plot we analyze the input vs. output SNR values of a double 
hop system assuming all four channels are Rician fading channels. 
The UHF uplink channel, the S-band downlink channel, the S-band 
uplink channel, and the UHF downlink channel are the channels 
which are considered. As one can see from the plot, various 
permutations of the 4 channels (specular to multipath power ratio 
of the fading channel is varied) are considered. 

These are just a few of the vast number of results obtainable 
with this simulation. Section 5 contains various sample runs 
illustrating some of the various possible results obtainable. 
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6.2. Results: Voice 


In addition to the analytical results, we were also able to 
use the simulation to process sampled voice and actually listen to 
it. In particular, the following tests were made on the voice: 

First, we put the sampled voice through the simulation with no 
fading, interference, or thermal noise present, or in other words, 
"pure FM" . The output voice sounded exactly the same as the 
input. We listened to the voice (still pure FM) with the 
pre-emphasis de-emphasis filters both off and on, and could here a 
slight improvement with the latter case. With the expander 
compander on or off, difference in voice quality was to slight to 
be heard . 


Next, with no other disturbances present, we listened to the 
voice as a function of the thermal noise power. Results could be 


summarized as follows: 

Carrier-to-Noise Speech Quality 

Power Ratio (dB) 


+40 
+ 35 
+30 
+ 25 
+ 20 
+ 16 
+ 14 
+12 
+ 10 
+08 
+06 
+03 
-02 


Excellent 
Excellent 
Excellent 
Excellent 
Excellent 
Very good 
Good 

Good (very light static) 

Fair (light crackling) 

Fair-Poor (crackling louder) 

Poor (crackling loud) 

Very Poor (Basically only crackling) 
No sign of voice (only loud crackling 
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As the above indicates, the quality of voice went from 
excellent just above threshold, to poor just below threshold, 
consistent with our analytical results as well as theoretical 
results . 

We then listened to the voice as a function of the specular- 
to-multipath power ratio (dB) of the fading channel (see sample 
run 5.1 for channel parameters! , with no other disturbances 
present (input carr ier-to-noise power ratio at +40dB) . Within the 
range of s^ecular-to-multipath power ratios from dB and better, 
output voice quality was excellent. Within the range from -10 dB 
through +5 dB, voice quality was good (tonal characteristics of 
the voice were altered somewhat). Below -12 dB, voice quality was 
fair (tonal characteristics altered), with some crackling present. 

The above examples are a few of the many possible tests on 
voice. Any of the 5 scenarios can be tested on the voice, with 
all possible parameter variations. (note that the above tests 
correspond to scenario 1). 
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7. Proposed Research for the Forthcoming Year 

Because of the present system's very modular structure, 
simulation of other modulation techniques over any one of the five 
scenarios can be done with relative ease. Since the simulation is 
designed in such a way that the inphase and quadrature phase 
components of the input signal are being passed through the 
components of the system in a parallel fashion, any type of 
modulation scheme which can be broken up into these components can 
be easily implemented » For example, the present simulator could 
be used to test the performance character istics of duo-binary FM 
or tamed FM, both which use Jess bandwidth than standard FM. 

At the same time, one could use the simulation to process 
voice using various modulation techniques to see which types 
produce the best subjective voice quality. In the report, we 
propose to construct a system to allow for simple testing of 
actual voice input. This system will allow for simple interfacing 
of voice input to the IBM system on which the actual simulation 
took place. We plan to build this device using a personal 
computer as the controller. 

Also because of the systems modular structure, it would be 
possible test other types of components or noise sources. 
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